Off-equilibrium dynamics in a singular diffusion model 
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We introduce a schematic non-linear diffusion model where 
density fluctuations induce a rich out of equilibrium dynam- 
ics. The properties of the model are studied by numerical 
simulations and analytically in a mean field approximation. 
At low temperatures and high densities we find a long off- 
equilibrium glassy region, where the system evolves out of an 
initially pinned state showing aging and a slow decay of the 
autocorrelation as an enhanced power law, along with strong 
spatial heterogeneities and violation of the fluctuation dissi- 
pation theorem. 
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As fluids are supercooled below the melting temper- 
ature their structural relaxation becomes very slow and 
may result in a glass transition characterized by a dra- 
matic increase of the viscosity. Under these conditions 
the dynamical evolution of glass-forming liquids shows 
a markedly out of equilibrium behaviorcl. In the main 
approach to the glassy dynamics, the Mode Coupling 
theoryH, the dynamical equations are solved by resum- 
ming a non trivial set of diagrams. This theory predicts 
an equilibrium relaxation time To which diverges as a 
power law .-at a dynamical transition. Today it's well 
established^ that the quoted theory applies in a region 
located well before the point of structural arrest, where 
the relaxation time is found experimentally to diverge 
according to the Voghel-Fulcher lawEl 



r ~ exp[v(p c - p) *] 



(1) 



Therefore we do not have information from the theoret- 
ical point of view on the kinetics close to the dynamical 
transition where standard theories do not apply. In this 
paper we introduce a schematic diffusion equation with 
a phcnomcnologically chosen mobility which reproduces 
the equilibrium relaxation time (Q) observed in this re- 
gion. We then study the consequences of such an as- 
sumption on the out of equilibrium dynamics. Due to 
its relative simplicity the model is amenable of both nu- 
merical and analytical investigations allowing a complete 
description of its features. Despite the simplicity of the 
equilibrium properties, the out of equilibrium behavior 
is complex and in qualitative agreement with the known 
properties of real systems. Our analytical results could 
in principle be tested quantitatively in experimental data 
and in molecular dynamics simulations. 

Glassy systems are usually schematized as composed of 
particles rattling inside "cages" of typical size a formed 



by the neighbors. Although diffusion is unimpeded in- 
side the cells, motion over larger distances is strongly 
suppressed at high molecular densities because a global 
rearrangement of many particles is required. This glassy 
behavior is known to, be reproduced also in the simple 
case of hard spheres!. For suitably coarse graining of 
space and time scales, the dynamics of supercooled liq- 
uids is Brownian in its microscopic origin. Therefore we 
consider a diffusion equation for the variable p(f, t) which 
represents the coarse grained particle density over dis- 
tances of order a 



dp(?,t) 
dt 



M{p)V 



SF{p} 



5p 

In Eq. (|) F[p] = J dr[p\np + (1 - a) ln(l - p)\ is the 
entropy of the lattice gas-like modelEl, which we con- 
sider for sake of simplicity. More realistic forms of F[p], 
which take into account the attractive interactions be- 
tween particles can also be considered; here we find, how- 
ever, that the qualitative features of the model do not 
changea. M(p) is a mobility, specified below, which is 
supposed to capture the main features of the constrained 
cooperative dynamics of a dense fluid. r)(r,t) is a gaus- 
sianly distributed random field, representing the thermal 
noise, whose expectations are givenQ by < r](f,t) >= 0, 
< r)(r,t)r]{r* ,t') >= -2TV{M(p)V [6{f - f)6(t - 1')]}, 
where < • ■ ■ > means ensemble averages and T is the 
temperature in units of the Boltzmann constant ks- 

Now we specify M(p). As shown below, from Eq. (Q) 
the characteristic equilibrium relaxation time To behaves 
as M ~P being the average density, for low T. Then, 
in order to reproduce the behavior (fy) of To, we assume 
a local mobility of the formEl 

-i] 



7?(f, t) 



(2) 



M(p) 



— e v[p(f,t)- 



(3) 



where a rescaled density, so that p c — 1, has been con- 
sidered. Due to its generality, Eq. (||) is also suited 
for the description of different physical systems where 
a constrained cooperative dynamics is believed to play a 
fundamental role, such as granular materials!! We have 
studied Eq. (||) for a temperature quench, by simulations 
and in a mean field approach. We sketch our numerical 
results before entering a detailed mean field analysis. 

Eq. (^) has been simulated by a standard first or- 
der Euler discretization scheme on a 128x128 two- 
dimensional square lattice with periodic boundary con- 
ditions, starting from an uncorrelated high temperature 
initial state. The system is quenched to a very low tem- 
perature (results will be presented for T = 1CP 4 , but 
similar behaviors are found for different temperatures). 
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FIG. 1. R 2 (t) (a) and S 2 (t) (b) are plotted for a quench to 
T = 1(T 4 and densities p = 0.70, 0.80, 0.85, 0.90, 0.93, 0.95, 
0.96 (from left to right). The main figures show the outcome 
of the numerical simulation, the insets refer to the mean field 
approximation. The location of the three regimes described 
in the text are outlined by the numbers 1, 2 and 3. 

The decay of the average density fluctuations S 2 (t) = 
y/ < (p — p) 2 >, is plotted in Fig. l b . For low densities 
p a normal liquid region is observed; S 2 (t) very quickly 
decays to the constant value characteristic of the equi- 
librium state. By raising the density, one enters a differ- 
ent region. Here the vanishing of M slows the dynamics: 
this produces a long-lasting off-equilibrium glassy behav- 
ior before the equilibrium state is reached. The behavior 
of S 2 (t) shows that for large densities the dynamics can 
be divided in three regimes. Initially, for t smaller than 
a characteristic time r p , S 2 (t) remains constant. This is 
the first regime. An analysis of the system configuration 
p(r, t) in this time domain does not show any appreciable 
evolution: the system is pinned in the initial configura- 
tion. Then, for t > t p a second regime is entered and 
less dense regions equilibrates whereas high density zones 
are still practically frozen. This is characterized by the 
decrease of S 2 (t). In this regime, that will be referred 



to as slow evolution, one observes pronounced correlated 
spatial heterogeneities in the system. This spatial pat- 
tern is outlined by a slow decay, as a function of k, of 
the structure factor C(k,t) =< p(k,t)p(—k,t) >, that is 
consistent with a stretched exponential fit (see Fig. ||) 
C(k,t) ~ Cexp — [l(t)k] 2 ^ with p ~ 1/6 (at variance with 
the Gaussian decay of standard-diffusion), similarly to 
some experimental observational!! Finally the system 
enters the equilibrium state characterized by a constant 
value of S 2 (t) . This whole pattern is reflected by the be- 
havior of the particle mean square displacement R 2 (t), 
shown in Fig.l a , calculated through 



R 2 (t) = [ D(t')dt' 
Jo 



(4) 



where D(t) =< M(p) > is the average mobility^. Fol- 
low densities R 2 (t) ~ t in the whole time domain, as ex- 
pected for simple diffusion. As the density is increased 
toward the limiting value p = 1 three regimes can again 
be distinguished. After an initial linear increase (regime 
1) a progressively more pronounced inflection is observed 
in an intermediate time domain (regime 2) whose dura- 
tion is enhanced as p is increased. The same pattern is 
observed in both spin-glass like lattice,gas and Lennard- 
Jones molecular dynamics simulation!!! Asymptotically, 
in equilibrium, R 2 {t) ~ t, as for simple diffusion. 

In order to get analytical results we now introduce an 
approximation on Eq. (^) by first expanding, the loga- 
rithm on the r.h.s. of Eq. (||) to lowest orderliil and then 
by replacing the mobility M(p) with the effective diffu- 
sivity D(p) =< M(p) >. Since average quantities do not 
depend on the position, due to space homogeneity, one 
has D(p) — D{t). Eq. (||) then becomes 



dp(r,t) 
dt 



D(t)W 2 p(r,t) + f](f,t) 



(5) 



where the rescaling t — > t/[p(l — p)\ and T — > T = 
Tp(l —~p), has been performed, and < r](f,t)r](r',t') >— 
-2TD(t)V 2 [<5(r-r')d(t -£')]■ Transforming Eq. (§) 
into reciprocal space, one obtains the following for- 
mal solution for the two time correlator C(k,t',t) =< 
p(k,t')p{-k,t) >, (t>f) 



C(k,t',t)=e-l R2 ^ +R2 W k '- 
|c*(fc,0,0) +f 



= 2_R 2 (t')fc 2 _ 



(6) 



The whole problem is now reduced to the knowledge of 
R(t) which must be calculated self-consistently enforc- 
ing Eq. (||), where D(t) is given by D(t) =< M(p) >= 
Jq M(p)P(p)dp. Here P{p) is the probability distribu^ 
tion of the density field that, for Eq. @ can be show: 
to be Gaussian at all times. Then we have 
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D(t) = [2icS 2 (t)] 



-1/2 



M(p)e 



-(p-p) 2 /[2S 2 (t)] 



dp (7) 



2 



The quantity S(t) can be computed as S 2 (t) = 
(27r)~ d J k<A C(k, t, t)dk where A is a momentum cutoff 
of order a -1 . From Eq. (^]), S(t) is a function of R(t): 

S 2 (t) = h[S 2 (0) - qf]R- d (t)<i> d [V2AR{t)} + qf (8) 

where $ d [x] = f*y d ~ 1 exp(~y 2 )dy, 

q = (E d /(i)[A/(27r)] d , h = [d/{AV2) d ] and E d is the 
surface of the d-dimensional unitary hypersphere. No- 
tice that the asymptotic value of the density fluctuations 
S 2 (oo) = qT~p(l — ~p) vanishes at the point of structural 
arrest. 
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FIG. 2. The slow decay, as a function of k, of C(k,t) in 
the second regime, is plotted here against fc 2/J , with fx = 1/6, 
for increasing times. C(k, t) is consistent with a non Gaus- 
sian fit: C(k,t) ~ Ce" (i(t)fc)2M . The inset shows the de- 
cay of C(k,t',t) as a function of t — t in mean field for a 
quench to T = 10 -4 with p = 0.95. Different t' are shown 
(f = 1, 10 4 , 10 5 , 10 6 , 10 7 , 10 9 , 10 11 , 10 13 , from left to right; the 
last two curves collapse because the stationary state is entered 
and time translational invariance is obeyed). 

Eqs. (||, 0, ||) are a closed set of equations that can 
be studied analytically. In the present approximation 
the non linearity of M(p) is accounted for by a self- 
consistency prescription for the calculation of D(t). Sim- 
ilar approximation techniques are well developed and 
widely used in severaLJjield of statistical physicsEi pro- 
ducing reliable resultalil. From Eq. (^) the normalized 
correlator C(k, f , t) = C(k, t', t)/C(k, t', t') is given by 



<7(fe,^i)=e-[* 2 W-*W 2 



(9) 



showing that a scaling form C(k,t',t) = S[4>(t)/(f>(t')], 
with (j)(t) = exp{— R 2 (t)k 2 } 7 is obeyed as-suggested by a 
scaling approach to dynamical processestM 

An important issue to understand the off-equilibrium 
dynamics is the relation between the response function 

to a small perturbing field h^, t) = J* dr 5 ^. k ^j > , 

and the correlation function in the unperturbed situ- 
ation, C{k,t',t). In equilibrium systems, where the 



fluctuation-dissipation theorem holds, the quantity X = 
Tdx%(t' ,t)/dC(k,t' ,t) is equal to one. Out of equilib- 
rium this relation is violated, ajid, generally, A is a func- 
tion of t' and t: X = X(t',t)l3. In the present mean- 
field approximation, one may interestingly show that 
the generalized "fluctuation-dissipation" ratio (FDR) X 
is a function of the sole t': X(t') = {[<?(£, 0,0) - 
f] exp (-2k 2 R 2 (t')) + f}- 1 . Only if t' -> oo, the usual 
version of the FDR with X = 1 is recovered (notice that 
A < 1). 

In the following we will report the main results of. the 
mean field analysis referring to a longer publicationo for 
all the details. From the solution of the model one sees 
that if the density p is small or the temperature T is high 
one immediately enters the asymptotic stationary state 
that will be described later on. For high densities and low 
temperatures, on the other hand, the evolution remains 
markedly far from equilibrium for a long period and three 
dynamical regimes are found corresponding to different 
behaviors of R(t) (see Figs. 1), as discussed below. 

Regime 1 - Pinning: For short times, such that 
R{t) < A~\ as shown by Eq. (§), C(k,t',t) is essen- 
tially constant since \k\ < A and the system looks pinned. 
In this regime we also have <&d[v2Ai2(t)] ~ R(t) d , con- 
sequently S(t) ~ 5(0) and D(t) ~ D(0). Therefore 
R 2 (t) ~ D(0)t (see inset Fig. l a ). The duration of this 
regime is t p ~ A _2 Z3 _1 (0). Physically r p corresponds to 
the time the particle spends inside its cell (cage). 

Regime 2 - Slow evolution: Pinning lasts up to t p . 
For t > T p , we have R(t) > A, thus particles diffuse out of 
the cages and the evolution starts. For sufficiently long 
times, computing R(t) through Eqs. (|[ 0, |)) one finds 
R 2 (t) ~ b(lnt) s (see inset Fig. l a ), where b = const, and 
8 = 6/d. Eq. (||) implies that, for fixed t' the correlator 
decays as an enhanced power law (see inset Fig. ||) 

C(k,t',t) = exp{i? 2 (i')A: 2 }exp{-6[ln(i)] <5 fc 2 } (10) 

and S(t) ~ (lni) -3 / 2 (see inset Fig. l b ). A logarithmic 
relaxation of the density fluctuations is also observed in 
Molecular Dynamics—simulations of out of equilibrium 
liquid glass formersEZl. When also t' > r p , one has 
C(k,t',t) = exp {-b([Mt)] s - [\n(t')] 5 )k 2 }. The char- 
acteristic duration time, r e , of the slow evolution regime 
can be estimatedQ, at low T, to be r e ~ M^i^p). 

Regime 3 - Stationary state: For long times 
t > r e a simple diffusive behavior is obtained because 
D(t) always attains asymptotically a constant value 
D(oo). This imply R 2 (t) ~ D(oo)t, as can be seen in 
Fig. l a , so that the normalized correlator exhibits the 
usual exponential decay as a function of t: C(k,t',t) = 
e -D(oo)tfc 2 e {K 2 (t')fc 2 } ( see inset Fig |)_ when fJ > Tej 
we have C(k,t',t) = e -D{^)[t-t']k 2 and time trans i a . 
tional invariance is obeyed. In the small temperature 
limit the densitii-.fluctuations 5(oo) can be approxima- 
tively neglectedt3 and D(oo) ~ M(p). This leads to 
t = M _1 (p) which gives Eq. (pi), as previously stated. 
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So far we have studied the out off-equilibrium evolu- 
tion of a system governed by Eq. (§) in the presence of 
a vanishing mobility for which Eq. (fil) holds in equilib- 
rium. We also want to consider the case in which the 
mobility vanishes as a power law, as is found for instance 
in the Mode-Coupling Theory of supercooled liquidsH: 
M(p) = [1 — p(f, i)] 7 . This relation implies an algebraic 
divergence of To: To = M(p) , as experimentally found 
in supercooled liquids in the temperatures or densities 
regions far away the ideal glassy transition. The dif- 
ferent form of the mobility, as stated before, does not 
changes the global picture described so far with three 
different regimes. In the first and third regimes R 2 (t) is 
linear in t as before, while in the second one we find, for 
7 > I, an anomalous diffusion R(t) ~ wt 13 , with w n — 
const, and j3 — 4/ ('yd + 4). In this regime we also findta a 
stretched exponential decay of the normalized correlator 
C(k, t', t) ~ exp {R 2 {t')k 2 } exp {-wt^k 2 }. 

In this paper we have introduced a phenomenological 
equation for off equilibrium glassy dynamics. The only 
ingredients of the model are the diffusive behavior and 
the request of a Vogel-Fulcher (or algebraic) divergence 
of To, obtained by assuming a mobility as in Eq.(^). With 
these sole ingredients the out of equilibrium evolution of 
the model is observed to be highly non trivial, even in 
the mean field approximation which we have studied in 
details. Consistently with mean field theory, also the nu- 
merical integration of the full model shows the existence 
of a gradual crossover from a normal liquid to a glassy 
behavior by raising the density. Some properties that are 
observed in systems close to the glassy transition, such 
as the existence of strong spatial heterogeneities, anoma- 
lous diffusion, slow decay and aging of density autocor- 
relations, are exhibited by the model. These predictions, 
as long as the non trivial fluctuation dissipation ratio 
X(t'), are all amenable of experimental check. In mean 
field this whole richness is observed in the preasymptotic 
off equilibrium dynamics (which, however, may be expo- 
nentially long), whereas the asymptotic equilibrium evo- 
lution is trivial. This is an important difference with real 
glassy systems, where a non trivial decay of C(fc,i',i) is 
also observed in equilibrium. However in mean field a 
non exponential decay of C(k,t',t) can be ruled out on 
general grounds due to the Doob's theorerr£3. Further 
studies are in progress in order to characterize the com- 
plicate fluctuations occurring in the equilibrium state. 
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